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Abstract 

The paper deals with classical problem for cracks dislocated in a cer- 
tain very specific porous elastic material, described by a Cowin-Nunziato 
model. We propose a method based upon a reducing of stress concentra- 
tion problem for cracks to some integral equations. By applying Fourier 
integral transforms the problem is reduced to some integral equations. 
For the plane-strain problem we operate with a direct numerical treat- 
ment of a hypersingular integral equation. In the axially symmetric case, 
for the penny-shaped crack, the problem is reduced to a regular Fred- 
holm integral equation of the second kind. In the both cases we study 
stress-concentration factor, and investigate its behavior versus porosity 
of the material. More in particular the stress concentration factor in the 
medium with voids is always higher, under the same conditions, than in 
the classical elastic medium made of material of the skeleton. Further, 
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as can be seen, the influence of the porosity becomes more significant for 
larger cracks; that is also quite natural from a physical point of view. 

1 Introduction 

There are a number of theories about mechanical properties of porous materi- 
als. One of them is a Biot consolidation theory of fluid-saturated porous solids 
[1]. Typically, these theories reduce to classical elasticity when the pore fluid 
is absent. This is why Cowin and Nunziato proposed a new theory to describe 
properties of homogeneous elastic materials with voids free of fluid [2]. This 
theory is a special case of microstretch elasticity of Eringen, when micropo- 
lar effects are discarded. Eringen's microstretch theory is more appropriate 
for geological materials like rocks, soils since, this theory takes into account 
the intrinsic rotations and stretch of materials [5], [6]. These effects are not 
considered in Cowin-Nunziato theory. However, the last theory is mathemati- 
cally simpler to deal with, and it takes into account voids in porous materials 
without any other phase, like liquid and gas. 

A general theory of such materials is currently well-developed by many 
authors (see, for example [3], [4]), but too few concrete problems are solved, to 
allow one to estimate the practical merits of this model. 

Generally, this theory is founded on the balance of energy, where presence 
of the pores involves additional degree of freedom, namely, the fraction of el- 
ementary volume. As a consequence, the bulk mass density is given by the 
product of two fields, the void volume fraction and the mass density of the 
matrix (elastic) material. 

An exact explicit solution is well known in the problem for a crack (both 
linear and penny-shaped) dislocated in the classical linear elastic space [9] . If a 
normal load is applied to the faces of the crack, then the shape of the faces near 
the edge of the crack under this stress can be represented explicitly as a root- 
square function. This permits analytical calculation of the stress concentration 
coefficient in the classical case. Obviously, stress concentration analysis is also 
very important in engineering practice for porous materials. Therefore, the 
main goal of the present work is to construct a strict solution of the static 
crack problem for the line (plane-strain problem) and penny shaped (axially 
symmetric problem) cracks dislocated in the linear elastic granular (porous) 
space. 
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We first give a short summary of basic equations, then demonstrate appli- 
cation of the Fourier transform that allows us to reduce the problem to some 
integral equations and to construct a direct numerical collocation technique to 
solve this equation. In the plane-strain problem we treat numerically a certain 
hypersingular integral equation, and a special kind of the collocation technique 
can be also applied to such equations. Finally, we demonstrate in figures how 
the calculated stress concentration coefficient at the crack edge depends upon 
variation of some physical and geometric parameters. 



2 Governing Equations for Elastic Media with 
Voids: Plane-Strain Problem 



We refer the deformation of the continuum to a fixed system of rectangular 
Cartesian coordinates Oxyz. Let us consider elastic material with voids which 
possesses a reference configuration with a constant volume fraction u . The con- 
sidered theory asserts that the constant mass density g has the decomposition 
[2] g = 7*/, where 7 is the density of the matrix material, and v (0 < v < 1) 
is the volume fraction field. 

Let (0 = v — vq) be the change in volume fraction from the reference 
one. Then the linear theory of homogeneous and isotropic elastic material with 
voids is described by the following system of partial differential equations [2,8] 



where \i and A are classical elastic constants; a, (3 and £ - some constants 
related to porosity of the medium. Besides, u denotes the displacement vector. 
Obviously, if (3 — 0, then the elastic and the "porosity" fields are independent. 
Thus, in the case (3 = the stress-strain state is insensitive to the function 0. 
The components of the stress tensor are defined, in terms of the functions u 
and 0, by the following relations (<% is the Kronecker's delta) 



ji A u + (A + ji) grad div u + (3 grad = 



(2.1) 
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Let us formulate the plane-strain boundary value problem for the porous 
(granular) medium. In this case u = {u x (x,y), u y (x,y), 0}, and the basic 
system (2.1) can be rewritten as follows 

d 2 u x . o d 2 u x 



+ (1 - c 2 ) 



1-c 2 




d 2 u y | H d</> = Q 
dxdy dx 



d 2 u x d<fi 
dxdy dy 







(2.3) 



'du x du 



with 



cr. 



dUrr 



A + 2/i dx 



+ (1 



\ dx dy J 

2c 2 ) ^ + H<P, 
dy 



a 



yy 



A + 2/i 



= (1 - 2c 2 ) 



du 1 + du 1 + H ^ 
dx dy 



(2.4) 



-> xy 



du T du,, 



dy dx 

The positive physical parameters (2.3)-(2.4) are introduced as follows 

H TT 2 _ a 2 _ a 



H 



I 2 



(2.5) 



X + 2/i A + 2/i 1 (3 

where the first two numbers c , H are dimensionless and the quantities l\ , l 2 
have the dimension of length. 

Let us consider a thin crack of the length 2a with plane faces, dislocated over 
the segment —a < x < a along the x-axis. Let the plane-strain deformation of 
this crack be caused by a constant stress a yy = a applied at infinity (i.e. at 
y = ±oo). Then, due to linearity of the problem, it is easily seen that both 
the shape of the crack's faces and the stress concentration at the crack's edges 
are the same as in the problem with a solution decaying at infinity and the 
following boundary conditions corresponding to the case when the normal load 
—<t is symmetrically applied to the faces of the crack and there is no load at 
infinity. 

For the last problem the boundary conditions over the line y = are 



Txy = , = (\x\ < OO) , 



CTyy = —00 (M < 0>) , Uy = (\x\ > d) . (2.6) 
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It is proved (cf.[7,10]) that the boundary condition for the function (2.6) 
directly follows from a balance principle. 

Let us apply the Fourier transform along the x-axis to relations (2. 3), (2. 4) 
and the boundary conditions (2.6). Then Eqs.(2.3) reduce to the following 
system of ordinary differential equations, with respect to images of the func- 
tions u x , u y , (p. Therefore, these become functions of the variable y only, with 
the Fourier parameter s instead of the first variable x (all Fourier images are 
denoted by respective capital letters, for all physical quantities): 



f c 2 U" -s 2 U x + (l- c 2 )(-is) Uy-isH§ = 

(i - c 2 ) (-i S ) u' x + u; - c 2 s 2 U y + H $' = 

isU x -U' y + ll<$>" - + * = 0. 



(2.7) 



Let us introduce the new unknown function g(x) as 

( g(x) , \x\ < a 
u y (x,0) = { , (2.8) 



then the boundary conditions (2.6) in Fourier images become 

U' x -isU y = 0, $' = 0, U y = G(s) , y = , (2.9a) 

where 

a 

G(s) = J <7(0e X p(^K- (2-96) 

—a 

General solution of the system (2.7) is constructed in accordance with a 
classical theory of ordinary differential equations (cf. [10]). Due to natural sym- 
metry of the problem with respect to x-axis, we give here only solution for 
y > 0, which has the following form (here all physical quantities of dimension 
of length, including components of the displacement vector, are written in a 
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dimensionless form being related to l 2 ) 
( isH \ 



( U x \ 

Uy 

V $ J 



D 1 



1-N 

Hq(s) 
1-N 

1 



+ D 2 



( i sign(s) \ 
1 




e -\s\y + 



/ 



1-N + c 2 . . M \ 

— — + iy s 

-is){N — 1 + c 2 ) 

y 

2Nc 2 



(2.10) 



V 



#(1 -iV-c 2 ) 



, q(s) = Vs 2 + 1- N, 



J 



where N = (l\/l\)H , < N < 1 is the so-called "coupling number". 
Obviously, the unknown constants D±, D 2 , D 3 should be defined from boundary 
conditions (2.9a). These yield the 3x3 algebraic system 



f SH " {8) D 1 + sD 2+ (1 - N)sig 1 s) D :i = 



1-N 



q(s) D 1 + 



N-l + c 2 



2Nc 2 I s| 



H(l -N -c 2 ) 



(2.11) 



Hq{s) 



D ± + D 2 



I 1-N 
whose principal determinant is 



G(s), 



(1 — N)q(s)sign(s) 
N-l + c 2 



(2.12a) 



The particular ones (related to the three unknowns when one applies the 
Cramer's rule) are: 

, 2Nc 2 s\s\ 
Ai = -G(s) 



G(s) 



H(N - 1 + c 2 ) ' 
q{s)sign(s)[(l - N) 2 + 2Nc 2 s 2 } 
(1 -N)(N- 1 + c 2 ) 
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(2.126) 
(2.12c) 



A 3 = -G(s)q(s)s . 



(2.12d) 



Thus, the Fourier image P(s) of the normal stress a yy over the line y = as 
follows 



P(s) 



Hs 2 



2fi 1 — N 

G(s)\s\ 



D 1 + \s\D 2 + 



D» = 



'l-Nyq(s) 



N-l + c 2 
2Nc 2 s 2 (q - \s\) + (1 - N)(l - N - c 2 )q 



(2.13) 



By applying Fourier inversion to relation (2.13), with the use of a convolu- 
tion theorem, one can express the normal stress on the line of symmetry y = 
in terms of its opening: 



(i-AQ ; 

2/i 



6 

a yy (x, 0) = J g(0K(x - Z)d£ , (& = f ) > ( 2 -14) 

—6 



where 



1 00 I I 

tf(a;) = — / ^-[2Nc 2 s 2 (q-\s\) + (l-N)(l-N-c 2 )q\e- isx ds 
2tt J q{s) 



1 °° 

— / L(s) cos(sa;)ds, q = q(s) = V s 2 + 1 — N , 

7T 7 



(2.15) 



L( S ) 



g(s) 



2Nc 2 s 2 (q -s) + (l-N)(l-N- c 2 )q 



Finally, the boundary condition a yy (x,0) = —a , (\x\ < a) from (2.6) 
yields the main integral equation for the considered plane-strain problem with 
the convolution kernel 



V 

jg(Z)K(x-Z)dt = -(l-N) 2 ^-, \x\<b 



(2.16) 
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3 Reduction to Integral Equation in the Axi- 
ally Symmetric Case 

In this case the general equations (2.1) reduce, instead of Eqs.(2.3), to the 
following system of partial differential equations for the components of the 
displacement vector in cylindrical coordinates u = {u r (x,y), 0, u z (x,y)} 

/ 8 2 U r + 1 dUr_ _ Ur\ + ^ d 2 U r + _ ^ d 2 U z | ^ 8<j) _ Q 

y <9r 2 r cV r 2 y <9z 2 <9r<9,2 cV 

( 1 - c 2 ) — (— + — \ + c 2 ( + ldu A + + = o (3 1) 

\ dr r J \ dr 2 r dr J dz 2 dz 

i I ,g r 2 r q t g z 2 \ q \ dr r dz J 



with 



+ (1 - 2c 2 ) + — ± )+ H<j> 



A + 2/x dr \ r dz 

[l-2c 2 ) (^L + ^+^l + Hcf) , (3.2) 



A + 2/j, \ dr r J dz 



a rz du r du z 



fj, dz dr 

If a round thin plane crack is placed at the center of the cylindrical coordi- 
nates: z — , r < a and its faces is under a normal constant load — <t , then 
the boundary conditions over the plane of symmetry z = are 

d(j) 

Orz = , — = (r < oo) , o zz = -o- (r < a) , w z = (r > a) . (3.3) 
Let us apply the Hankel transform along the r-variable to relations (3.1), (3. 2): 

oo oo 

U r (s,z) = J u r (r, z) J\{sr) rdr , u r (r,z) = J U r (s, z) J\{sr) sds , (3.4a) 



oo oo 

V.(., z) = J Mr, «) Msr) rdr , «,(r, ,) = / (/,(«, ,) J («r) srfs , (346) 
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0( S , z) = J S(r, z) J (sr) rdr , $(r, s) = / *) J («r) . (3.4c) 
o o 

If one applies the Hankel transform to equations (3.1), then one comes to 
the system of ordinary differential equations, with respect to images of the 
functions u r , u z , 0. Therefore, these become functions of the variable z only, 
with some parameter s: 



< c 2 U" -s 2 U r -(l- c 2 ) sU' z - = 

(1 - J)sU' r + U'J - c 2 s 2 U z + H <!>' = 

- sUr -u> + il&>-( l l + il s A $ = 0) 



(3.5) 



where all derivatives are applied with respect to the variable z. 

Solution of Eqs.(3.5), together with Eqs.(3.4) and boundary conditions (3.3) 
gives for the Hankel image of the normal stress a zz the following expression in 
terms of the image of crack's faces opening 

P(s) G(s)s 



[1-N)*q(s) 



2Nc 2 s 2 (q - s) + (1 - N)(l -N- c 2 )q 



where 



and 



u z (r,0) 



g(r) 
, 



r <b 
r > b 



(3.6) 



(3.7) 



G(s) = J g(p) J (sp)pdp , P(s) = J a zz (p,0) J (sp)pdp . (3.8) 
o o 

If one applies inverse Hankel transform to relation (3.8) then one arrives at 
the main integral equation for the round crack problem: 

6 

J g(p)pK(r,p)dp = -(l-N) 2 — , 0<r<6 (b = f ) (3.9a) 



// 



with the kernel 



K(r, p) = J L(s) J (rs)J (ps) s 2 ds. 



(3.96) 



9 



4 Properties of Integral Equations and Numer- 
ical Implementation for the Plane Problem 



Let us start from the evident estimate 



q(s) ~ s 



1 - N 1 , 

1 + — + V- 



at s -> +oo , (4.1) 



so the asymptotic behavior of the symbolic function of the kernel at infinity is 
given as follows 

L(s) ~ (l-iV) 2 (l-c 2 )s + oQ , s^oo. (4.2) 

It is obvious that the leading asymptotic term (4.2) leads to a principal degree 
of the kernel's (2.15) singularity and one can see that the latter is defined by 
the integral 

oo 

/ scos(sx)ds = , (4.3) 

J X 


i.e. the kernel is hypersingular when x — > [9]. It can be shown that the 
remaining term in asymptotic estimate (4.2) possesses a logarithmic (i.e. weak) 
singularity at x — > 0. 

However, the kernel (2.15) admits explicit representation that permits direct 
estimate of its singular properties. To obtain such a representation, we calculate 
the following integrals 

oo l0 oo 

/dp 6 
s 3 cos(sx)ds = — -j-^ J scos(sx)ds = — . (4.4) 

o o 

Then we consider the following table integral (here K is a McDonald's function 
of the order 0) 

oo 

I cos^j^ = Kq{VY ^N\x\) , (4.5) 
and apply, step by step, sequential derivatives to Eq.(4.5), as follows (assuming 
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x > 0) 



DC 

I 



s 2 cos(sx) , d 2 
-as = — 



dx 2 



d 



d 



= Vl - N— KAVl -Nx) — (1 — N) — , 



K^VT^Nx) = (4-6) 



1 - N) K (VT^N x) + 



x 



x 



s 4 cos(sx)ds 



r s cc 

J (T^ 



(l-N)q(s) dx 2 



3 1-N\K 1 (VT^N 



d2 f K (VT^N X ) + *(^*) 



(4.7) 



+ 



+ 1-7V+ — X (v / T 3 ^V 



vX" a; / y/l — N \ x* 
where we have used the derivatives of the McDonald's functions of the zero's 
and the first order: 
dK (z) 



(4.8) 



dz dz z 

Therefore, the full kernel for the plane-strain crack problem is explicitly 
represented as follows 



K(x) = - I 2Nc 2 



x 4 



l-N 









\x\ 



K^Vl- N\x\)- 



(1 _ N ) ( i _ N + ) Koiv^N 



x/ 



x 



1 - N)(l - N - c 2 ) 



(4.9) 



x z 



^From the last representation it is obvious that the function K(x) (4.9) is 
differentiable as many times as wanted outside a small vicinity of the origin 
x — 0. Let us estimate its behavior when x — > 0. For this aim one can use the 
following asymptotic formulas 



K Q (z) ~ -In^-j - 1 + 0(z 2 \nz) , z - +0 , 
~ i + £ln(|) + ( 7 -i)| + 0(^lnz), ^ +0 , 



(4.10) 
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where 7 = 0.577216 is the Euler's constant [10]. Thus behavior of function 
(4.9) for the small argument is 



K(x) 



(l-N) 2 (l-c 2 ) 



7T X" 



+ 0(ln|x|) , x -> , 



(4.11) 



so the kernel is hypersingular indeed. 

The numerical method that we apply to solve the hypersingular integral 
equation has been proposed in our previous paper [9]. It is based on extraction 
of a characteristic hypersingular part of the kernel. Let us represent the full 
equation as follows 



u 

J 



1 



(x-ty 



+ K* (x, t) 



g (t) dt = f'(x), x G (-6, b) , 



(4.12) 



where the regular part may admit a weak (i.e. integrable) singularity. If 
we represent it as 

dK 1 (x, t) 



(x, t) = 



dx 



(4.13) 



then a bounded solution can be constructed by applying inversion of the charac- 
teristic part, that reduces eq.(4.12) to a second-kind Fredholm integral equation 

6 

g (x) + J N, (x, t) g (t) dt = ft (x) , x G (-6, b) , (4.14) 



where 



TVi(rM) 



x" 



7T^ 



-6 



K l (r, t) dr 
\Jb 2 — t 2 (x — t) 

f(r)dr 

Vb 2 ^ 



T z (X — TJ 



(4.15a) 



(4.15a) 



Then we prove that, if / (x) G C\-b,b); K^x^) G C 1 [(-6,6) x (-6,6)], 
then for any x G (—6, 6) the difference between solution g (x) of the linear 
algebraic system 



E 

j=l L 



1 



X{ tj X% — 1 



+ hK* (xi,tj) 



g(t j ) = f'(x i ), i = l,...,n (4.16) 
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and the solution of eq.(4.14) tends to zero near the crack edge when h — > 
(i.e. n — > oo), where the mesh nodes X; L = — b + (i — 1/2) h, i = l,...,n and 
tj = — b + jh, j = 0, 1..., n are taken with the constant step h = 2b/n. 

It should be noted that this result automatically implies: if our numarically 
constructed solution of the system (4.16) tends to the exact one, so the con- 
structed solution tends to zero when approaching the crack edges, as it follows 
from (4.14)-(4.15). 

5 Penny-Shaped Crack: Reduction to a Fred- 
holm Integral Equation of the Second Kind 

Here we apply a special transformation rather typical for the problems with 
axial symmetry. 

Let us rewrite Eq.(3.9b) in the following form 

b oo 

Jg(p)pdpjL(s)J {rs)J (ps)sds = -(l-N) 2 ^, 0< r < b , (5.1) 





and introduce the new function as follows 



/ g(p)J (ps)pdp = J g(p)J (ps)pdp = J 7 (£)5^d£ . (5.2) 



The necessary and sufficient condition for Eq.(5.2) to be correct is that this 
must provide g(r) =0 for r > b. Let us control correctness of this statement 
by applying inverse Hankel transform to (5.2): 



g (r) = j J ( rs )ds J 7(0 sin(£s)d£ = 



5 , r < f 



r , r > f 



><%={ 



b 



, r > b 



(5.3) 



so (5.2) is correct. 
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Now Eq.(5.1) becomes 

b 

7(0 sin(£sR [ L(s) J (rs) ds = -(1 - N) 2 

fi 

u u 

and one may apply the operator 



U CXJ 

J 1 (Osm(^s)d^jL(s)J (rs)ds = -(l-N) 2 —, < r < b , (5.4) 



to the both sides of Eq.(5.4), that leads to the following integral equation with 
respect to the function 7(£) 



oc 



/ 7(£K / — sin(£s) sin{ts)ds = -(1 - N) 2 — , < t < b , (5.6) 



if the applied normal load cr Q is constant. When performing all above transfor- 
mations, we have used the table integrals 



7jb(r.) \7W^' r<i , r J^L = ^M. (5.7) 

Let us extend the uknown function 7(0 in Eq.(5.6) to a negative interval 
£ G (—b,0) as an odd function. Then the integral (5.6) can be represented as 
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follows (S(x) is the Dirac's delta) 



J 7(0* 



sin(£s) sin(is)ds 



cos[(£ - - cos[(C + t)s\ 



ds = 



1 ? 



7(0* 
L( S ) 



cos[(£ — t)s] ds 



(5.8) 



- (l-iV) 2 (l-c 2 ) 



cos[(£ — i)s] ds+ 



+|(1 - iV) 2 (l - c 2 ) 7 (t), since ^ cos(s0ds = ^(0 • 



Therefore, in this axially symmetric problem one arrives at a very regular 
Fredholm integral equation of the second kind with the convolution kernel: 



*i{t) + / i(0K*(t - 0* = - (1 , \t\<b, 



-b 



(5.9a) 



where 



00 

K*(x) = J 



L(s) 



(1- N) 2 (l-c 2 )s 



cos(a;s) ds 



(5.96) 



2Nc 2 s 2 (q - s) + (1 - N)(l -N- c 2 )q(s) 
(l-N) 2 (l-c 2 )q(s) 



- 1 



cos(sx)ds. 



Since the expression in the square brackets here is of the order 0(l/s 2 ) at 
s — > 00 (see (4.2)), it can be easily proved that the kernel K*(x) is regular. 
More precisely, it is differentiable over any finite interval, in particular K*(x) G 
C\-2b,2b). 
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6 Calculation of the Stress Concentration Co- 
efficient 



Here we are intersted in behavior of the normal stress (which is a yy for in-plane 
problem and a zz for axially symmetric problem) near the crack edge. 

Let us start from the problem with axial symmetry. We first notice that 
function 7(£) is regular on the interval £ G (—6, 6), as a solution of the regular 
Fredholm equation of the second kind (5.9). Then we operate with the relation 
(5.4), which obviously determines the normal stress a(r, 0) not only for < 
r < b, but also for r > b. It is easily seen that with r — > b + the leading 
asymptotic term is given by the leading term of L(s) at s — > oo: 



(J, 



(r,0) 



(1 - 
= -/7(0 



d_ 



o o 
1 



=5 , i <r 



^/r 2 - e 

0, £>r 



"/ 7 ^d£ / cos (£ s ) J oMds = 



/ 7(0^U= (6-1) 



7(0-7W 7W 



7(6) 



6 2 



6 2 



at 



6 + , 



so the stress concentration coefficient is 
k= hm |a - M)l 



r^b+O jxb{\ - C 2 ) 



6 2 = 



17(6)1 



(6.2) 



Some examples of calculation of this dimensionless coefficient for the penny- 
shaped crack are shown in Figures 1 and 2. 

In the plane-strain problem our approach is absolutely different, since so- 
lution of the integral equation (2.16) with the hypersingular kernel (4.9) is 
constructed numerically from the algebraic system (4.16), so that the applied 
numerical method does not operate explicitly with the root-square structure 
V6 2 — x 2 , vanishing at the edges of the crack. Therefore, we need to treat 
expression (2.14), which represents the normal stress a yy (x,0) also for x > 6, 
in a direct numerical way. 
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Since K(x) = -(1 - N) 2 (1 - c 2 )/(nx 2 ) + 0(ln \x\) , x -> 0, it is obvious 
that with x — > 6 + the leading asymptotic term of the considered normal 
stress is 

2MT^T = (i - ATO - ^ J gm{x - ^ ~ -«JTZ=W> (0) 



-6 



so here the stress concentration factor is 



A: 



hm %^liv^^2 = 2 Um 

x^b+O jlb{\ — C 2 ) 7T 



Vx 2 - 6 2 . (6.4) 



Figures 3 and 4 demonstrate numerical results on computation of this di- 
mensionless constant for some values of parameters N, c 2 , b = a/l 2 - 



7 Conclusions 

1. We have proposed a method based upon a reducing of stress concentration 
problem for cracks to some integral equations. In the plane-strain problem this 
is a hypersingular integral equation, which permits efficient direct numerical 
treatment. In the axially symmetric problem for penny-shaped crack, after 
some traditional transformations, we arrive at more regular Fredholm integral 
equation of the second kind. The both types of equations can generally be 
solved numerically. 

2. In some cases the developed equations admit exact analytical solution 
in explicit form. The first case is for N = 0, that is a classical linear elastic 
material. In this case the plane-strain problem's integral equation has only 
the characteristic component —(1 — iV) 2 (1 — c 2 )/(nx 2 ) and the regular part 

vanishes. As can be seen from Eqs.(4.15) here Ni = 0, so exact solution 
is given by (4. 14), (4. 15b), that is obviously coincides with a known classical 
solution. For the penny-shaped crack, the kernel K* also vanishes, as can be 
directly seen from Eq.(5.9b), so function jit) can be explicitly extracted from 
Eq.(5.9a), that finally also leads to a well-known classical solution for the round 
plane crack. 

Another simple limiting case is when parameter b = a/h is small, that 
means the crack size to be small when compared with the physical parameter 
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li (the latter being of dimension of length). Indeed, in the axially symmetric 
problem, if b — > 0, then integral operator in Eq.(5.9a) vanishes. The same 
property takes place for the plane linear crack, since integral of the regular 
part of the kernel in (4.12) vanishes. Thus, in the both cases the stress-strain 
state near the small crack in the porous space is like in an ideally elastic classical 
medium, that is quite natural from the physical point of view. 

3. If we investigate the influence of the porosity to the stress concentration 
factor shown in figures, we can discover very interesting properties. First of 
all, this factor in the medium with voids is always higher, under the same 
conditions, than in the classical elastic medium made of material of the skeleton. 
This can be explained with energetic arguments, since in porous media a stress, 
distributed near the crack edge only in the skeleton, can provide a balance of 
energy caused by applied normal force. So this requires more intensity of the 
internal stress to provide the balance. 

Further, as can be seen, influence of the porosity becomes more significant 
for larger cracks; this is also quite natural from a physical point of view. This 
results is very interesting conclusion. Let us imagine a relatively small crack, 
which by a sufficiently high stress concentration factor extends slightly in its 
length. Then, as clear from the figures, the larger crack is coupled with higher 
stress concentration factor that causes further crack extension. This is the real 
physical mechanics of cracks growth in porous media. 
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LEGENDS TO FIGURES 



Fig 1. Relative value of the stress concentration factor k with respect to 
ko in classical elastic medium versus coupling number, penny-shaped crack: 
c 2 = 0.2 . 



Fig 2. Relative value of the stress concentration factor k with respect to 
ko in classical elastic medium versus coupling number, penny-shaped crack: 
c 2 = 0.4 . 



Fig 3. Relative value of the stress concentration factor k with respect to ko in 
classical elastic medium versus coupling number, plane linear crack: c 2 = 0.2 . 



Fig 4. Relative value of the stress concentration factor k with respect to k in 
classical elastic medium versus coupling number, plane linear crack: c 2 = 0.4 . 
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Fig.2 




Fig.3 




Fig.4 



